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. Abstract 
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It is shown that constant galactic rotation curves require a logarith- 
I mic potential in both newtonian and relativistic theory. In newtonian 

04 ' theory the density vanishes asymptotically, but there are a variety of 

possibilities for perfect fluid einstein theory. 
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1 Introduction. 



Constant galactic rotation curves require a logarithmic potential in both 
newtonian and relativistic theory. The logarithmic potential contrasts with 
heuristic derivations of the Tully-Fischer [H [7| relationship, which require 
(j) ~ 1/r. As pointed out in |21j for a potential ~ ln(r) one would expect as 
well a metric with ~ ln(r) terms and hence first derivatives and Christof- 
fel symbol corresponding to the metric are of order ~ 1/r; thus both the 
derivatives and products of the Christoffel symbol are of order ~ 1/r^ so 
that linear approximations to the Riemann tensor and hence the field equa- 
tions are not necessarily consistent. This is why for a relativistic analysis 
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one has to start again with the properties of geodesies, this is done here in 
section [3l A problem with any model of galactic rotation is what happens 
asymptotically. The logarithmic potential necessary for constant rotation 
curves diverges as r ^ cxd. One way around a divergent potential seems to 
be to invoke a "halo" which is a region surrounding a galaxy where there 
is an approximately logarithmic newtonian potential. The properties of the 
"halo" have some sort of unspecified cut off so that the potential does go out 
to too large distances, quite how this works is not at all clear, after all binary 
and multiple galaxies exhibit unusual dynamics, as do stellar clusters, but 
presumably nearby galaxies have more gravitational influence than distant 
ones, so that at some distance a logarithmic potential must have a cut off. 
Another way of looking at what happens asymptotically for a logarithmic 
potential is to note that it is only the derivatives of the potential that have 
effects and (pr ~ ^/r does not diverge asymptotically. Yet another way is 
to introduce smooth step functions, see section [H "Halo" is a term that is 
usually restricted to an area of "dark matter" surrounding a galaxy; however 
if a model uses matter fields or different field equations to model galactic 
rotation one would again require a region that these fields model and it is 
simpler to again refer to this as a "halo" rather than use new terminology. 

In section m the analysis is generalized to the perfect fluid Einstein equa- 
tions where it is shown that the density is the same as in the newtonian 
analysis, but with additional terms appearing in the pressure and the den- 
sity of order v'^. The analysis requires four assumptions: l)a galactic halo 
can be modelled by a static spherically symmetric spacetime, 2)in the limit 
of no rotation halo spacetime is Minkowski spacetime, 3)the constraint on 
the first derivatives of the metric needed in order to produce constant ro- 
tation curves can be applied, 4)the perfect fluid Einstein equations; the 
density is the same as in the newtonian analysis, but with additional terms 
appearing in the pressure and the density of order v^. The assumption 1) 
is that a galactic halo can be modelled by a static spherically symmetric 
spacetime. That it is static is a good approximation as dynamical distur- 
bances in galaxies do not seem to be the cause of constant rotation curves. 
That it is spherically symmetric seems to be born out by observations [16]. 
Thus it is taken that a spherically symmetric halo causes constant galactic 
rotation, rather than the rotation of the luminous part of the galaxy. Using 
spherical symmetry as opposed to axial symmetry considerably simplifies 
calculations. There are two problems with the assumption 2) that in the 
limit of no rotation the spacetime is Minkowski spacetime. The first is that 
it excludes the possibility of other limiting structure, such as deSitter space- 
time. The second is that even non-rotating galaxies have non-flat spacetimes: 
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in other words this assumption is that to lowest order the rotating part of 
the description of a galaxies spacetime can be described independently of 
its non-rotating part. The assumption 3) is the constant velocity rotation 
constraint, this is a constraint on the first derivatives of the metric, it has 
been derived in |22j . and is derived again in the next section. For spherical 
symmetry it is an unusual constraint because it often involves logarithms 
rather than powers of the metric. The assumption 4) is the perfect fluid 
Einstein equations. For spherical symmetry the stress has three arbitrary 
stress functions {p,pr,pe)- The assumption that the stress is a perfect fluid 
reduces the three components to two {p,p) via the isotropic pressure equa- 
tion p = Pr = Pe- For a perfect fluid the density p and the pressure p 
are independent unless an equation of state is imposed. The reduction of 
the three arbitrary stress functions to two can be thought of as a second 
constraint on a first derivative of the metric. The fourth assumption has 
two aspects: firstly that Einstein's equations are the correct ones to model 
galaxies, and secondly that any "dark matter" has isotropic pressure so that 
Pr = Pe- In section U whether it is permissible to expand in the velocity of 
rotation Vc/c, where c is the velocity of light, rather than a radial param- 
eter r is also discussed. For a galaxy it is not clear what radial parameter 
could be used because of the problem of what happens as r — > oo, for the 
luminosity radial parameter R one has that the metric ~ \n{R) suggesting 
that it is not a good parameter to expand in. In both cases the expansions 
are independent of the mass of the galaxy m which does not enter explicitly 
into calculations here to lowest order. Taking galactic rotation to be in any 
way related to the Kerr solution is always misleading, as Kerr rotation is a 
short range effect and galactic rotation is a long range effect. 

A selection of recent literature on rotation curves can be divided into 
three types. The first are papers which use the MOND modification of 
gravity, see for example [21 [3] . The second are papers which invoke higher 
dimensions and/or branes, see for example [HI [lOl [TTl [151 [IS]- The third 
are papers which invoke new fields, or field equations, or potentials, see for 
example [l9l[Mll2ll[6l[l2l[l3l[ia[Il[5l[2ll[25l[26]. 

Conventions include: metric -+++ and field equations Gab = SnGTab 
Calculations were done using grtensor2/maple9 [17| . 

2 Newtonian theory. 

First a newtonian analysis is given, the presentation follows [22], but is 
shorter. In spherical coordinates newton's second law has two dynamical 
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components 

-{<Pr,M=(^r-re^^ir^)y (1) 

where (p is the newtonian potential and the subscript indicates which coor- 
dinate it is differentiated with. The 6 component integrates to give 

r^e = L, (2) 

where L is the angular momentum for each individual particle. The square 
of the velocity at any point is 

=r^ + {rOf. (3) 
For circular orbits r = so that ([3|) becomes 

= {re? = ^, (4) 
also for circular orbits newtons second law ([1]) gives 

cPr = r9' = ^. (5) 
Eliminating from these two equations and ^ gives 

r 

integrating gives the newtonian potential 

(^ = w2ln(r). (7) 

Poisson's equation is 



br = 



(6) 



AnGp = Oct> = -^{r'<Pr)^ = ^. (8) 

3 The geodesies of rotation curves. 

The constant rotation curve constraint can be derived as follows. The axi- 
symmetric line element can be taken as 

ds^ = gttdt"^ + 2gt^dtd(p + g^^dcp^ + grrdr"^ + geedO^- (9) 
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The geodesic lagrangian is 

where 2vC = — 1 for timehke geodesies. The euler equations are 



dpa _ dC 

dr dx°- 



fill 



For ^ 

dC dC 

Pr = ^ = 9rrr, 2— = g[^f + 2g[^i4> + g'^^^^ + g'^rV + g'ee^^. (12) 
where "' denotes d/dr. Thus the r component of the euler equation is 

^ = Qrrf + g'„r'^ + grrfirO = g'tft^ + 2g't^i(j) + g'^^^i^ + g'^^r^ + g'ggd'^, (13) 

with a similar equation for the 9 component. If required the pt and p^j, 
components can be expressed in terms of the energy E and momentum L 
for each geodesic, and then with [13] this ahows generahzations of the Binet 
equation; however for present purposes this is unnecessary as the problem 
much simphfies because constant velocity curves require 

e = r = 0. (14) 

Substituting into the lagrangain (jlOp and the r ^ 9 components of the euler 
equation (fT3l) gives 



2C = gut' + 2gt^t(t) + 5</,00^ (15) 
Q = g[,P + 2g[^i^ + g'^^4>\ 

respectively. There is a similar equation to the second of (|15p with'" denot- 
ing d/d9. The angular momentum is 



f^^^ = ^, (16) 



note that this corrects equation (27) of [22] where there is an r rather than a 
yjgcjxf, in the denominator of the last term. Substituting into the lagrangian 
equation of (jl5|) from which t can be calculated, this is not necessary for 
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present purposes. Substituting into (jlSp first for the first equality in ([16 
and tlien for the second equahty in (jl6p gives 



^ = g'u + 2g[^^ + g'^^^' = g[^ + 2-^gt^ + (17) 



For spherical symmetry with line element 

ds^ = - exp{2i')dt^ + exp{2fi)dr^ + exp(2V')fiS|, (18) 

where (iS2 = dO'^ + sin(0)^(i0^. The constraint (jlTh reduces to 

{exp(2z.)}' = 2t;2^'. (19) 

The requirement that for Vc = the metric is Minkowski fixes the constant 
of integration to be —1; integrating the metric becomes 

ds"^ = - (1 + 2v'^i;{r)) dt^ + exp{2fi{r))dr'^ + exp(2V'(r))(iS^. (20) 

At first sight this line element has two arbitrary function in it, however 
defining the luminosity distance r' 



r' = exp(^(r)), eMAr') ^ '^'''''''TP ' ^''^ 



exp(/i(invV'(ln(r')))) 

and dropping primes the line element becomes 

ds"^ = -(1 + 2v'^ ln{r))dt'^ + exp{2fi)dr'^ + r'^d^l, (22) 

leaving just one arbitrary function /i(r). This can be fixed by using field 
equations. 

4 The perfect fluid Einstein equations. 



For a spherically symmetric spacetime of the form (j22p there are three non- 
vanishing components of the Einstein tensor 

8^Gp = -G*i = ^(l-(rexp(-2^*))0, (23) 
STrGpr = = — ( -l + exp(-2^)^ ^ ^ ^ 



r2 V ' - (1 + 2't;2 ln(r)) 

-exp(-2/i) / 
r2(l + 2w2ln(r)) V (1 + 2^2 ln(r)) 



STTGpe = +G%= "Z^,^,^\A,, , , ,, +V(l + .g + 2.gln(r))^ 
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where p is the density, pr is the radial pressure, and is the angular pres- 
sure. These field equations are simple as they linear in only the first deriva- 
tives of /X and so should have many solution, for present purposes just the 
case of a perfect fluid is examined. The requirement that the stress be that 
of a perfect fluid is that the pressure is isotropic p = Pr = Pe, and this is 
sufficient to completely determine the metric. Subtracting pg from pr gives 
a first order differential equation in fj,, 

, l + vj + 2v'^.\n{r) , x n 1 + 2t;2 ln(r) 

(24) 

setting pq = pr and solving gives 

the requirement that the line element reduces to Minkowski when the veloc- 
ity vanishes gives = 0. / is the integral 

^ ^ j ^ {vl + l + 2vl\n{r)) ^^^^ 
= + (1 + ^ - + + 41n(r)2) ^ + {vl) . 

For A: = the pressure is 

SvrGp = -^ + 2[l + 2vl + 2vl ln(r)) (A; - /) = ^ + O [v^) , (27) 
and the density is 

MGp = (3 + vl + ln(r)) ^ g 2 _ ^) ^^g^ 

= +^-(9 + 4Mr))^ + OK«). 

The 0(1;^) term of the integral (j27p shows that the spacetime is Minkowski 
when Vc = k = for k ^ the (^^r component looks similar to deSitter 
spacetime however there is no corresponding term in gu and any such term 
would violate the the rotation constraint (jl7p . here usually k = 0. The 
expansion in Vc shows that to lowest order the density (j29p agrees with the 
newtonian expression however the expansion does not always converge, 
terms becoming larger and changing sign for each increase in v'^. 
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density 

pressure 

mass 



Figure 1: The density, pressure and mass function. 
The integral can be expressed in terms of the exponential integral 

^ = (l + ^) Ei (l' 1 + ^ + 21n(r)) , (29) 

and this can be used to numerically produce graphs of the properties of the 
spacetime. In figure one the density, 10^ x pressure and 400 x mass function 
are plotted for Vcjc = 0.05 (not v\ = 0.05) and k = 0. Asymptotically the 
density approaches 3, the pressure 0, and the mass function r. In the region 
Vc = 0.01 — 0.02 the machine overflows, which is unfortunate as Vc = 10~^ 
is a more realistic value. For k = 1 the density and mass function appear 
negative, the critical value for which positive density is restored is around 
k = 2.5x 10~^, there does not seem to be a critical value for the mass function 
which always goes negative for large enough r. For k = —1 the pressure is 
negative, the critical density for which positive pressure is restored is around 
k = -2.5 X 10~i°. 

The null convergence condition is RCC = Rabn°"n^ > 0, for a null 
vector p.95[9], the weak energy condition is TCC = TabV"'V'', for Va a 
timelike vector p.89[9], the dominant energy condition is the timelike con- 
vergence condition and T°'^Va is a non-spacelike vector p91[9], NSV is the 
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RCC, TCC and NSV for a rotation of 5% 
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Figure 2: The energy conditions. 



size of T'^'^Va. In figure two lO^xRCC and lO^xTCC and -lO^xNSV are 
plotted for Vc = 0.05 and A; = 0. For k = 1 the nuh convergence condition is 
violated, the critical value around which it seems to be restored is A: = 10~^. 
For k = —1 the weak energy condition is violated and this seems always to 
be the case for large enough r. The curvature invariant RiemSq is defined by 
RiemSq = Rabcd.R"''^'^'^ , with similar definitions for WeylSq and RicciSq. In 
figure three 10^ x RiemSq, 10^ x WeylSq, 10^ x RicciSq and lO^xRicciscalar 
(not squared) are plotted for Vc = 0.05 and k = 0. It does not seem to 
be possible numerically to determine whether the divergence happens at 
r = exp(— l/(2t;^)) ~ 10^^^ or r = or both. A surprising feature of figure 
three is that WeylSq is large compare to RicciSq, this has the interpreta- 
tion, see [9]p.85, that more of the curvature is due to gravity as opposed to 
matter. For k = ±1 the RicciSq is large compared to WeylSq, the critical 
value seems to be around k = ±10"^, but seems to highly dependent on r. 
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Curvature invariants for a rotation of 5% 
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Figure 3: The curvature invariants. 



5 Asymptotics and units. 

Constant velocity curves are only observed over a certain region, how short 
and long radial distances fit to this region is a problem. There seem to be 
three approaches to this: the first is to adjust things at the last moment and 
produce an onion model in which the constant velocity region as given has 
other spacetimes are fixed to it, see [23], the second is to adjust things in 
the middle by choosing contrived newtonian potentials, the third is to adjust 
things at the beginning by letting Vc Vc{r). 

Looking at the third approach first, for a newtonian model, equations 
([6]) and ([8|) become 

rct>r = 7;^(r), 4^Gp = ^ (rt;2(r))' . (30) 
For purposes of illustration consider the smoothed out step around r = ro 

Sir,ro) - ^ (1 - tanh(r - r.)) = ^^^(.(^ ! ^^)) ^ ^ ■ (31) 
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constant velocity galactic rotation curves with a step cutoff 
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velocity 
density 



Figure 4: A good newtonian potential. 



Apply the step ([^T]) to the velocity 

Vc{r,ro) = S{r,ro), 

then the potential is given by an integral, 
has to be evaluated numerically 



dr 



(32) 



with k = 1, The integral 



r(exp(2(r — Tq) + k) ' 



for the velocity (j32|) . (|30|) gives the density 



AirGp 



1 + (1 - 2r)exp(2(r - r, 
r2(l + exp(2(r — Tq)))^ 



(33) 



(34) 



For To = 4 the velocity (j32|) . the potential (|33|) divided by 3.5 and the density 
()34p divided by 100 are plotted in figure four, to avoid singularities at the 
origin the integral is taken from 0.1 rather than 0. By comparison with the 
integral / (ir/(exp(2(r — Tq) + 1) = ln(exp(2(r — ro))/(exp(2(r — ro)) + l))/2 
the potential converges to as r ^ oo, which is good for two reasons, firstly 
that it converges at all unlike ln(r), secondly that the limit is so that 
there is no trace of the potential at great distance from the galaxy; similarly 
comparing with the integral / dr/r = ln(r) the potential is seen to diverge 
for short distances. Surprisingly both the potential and the density show no 
unusual properties at the step Tq. 

For the second approach apply the step [3l] to the potential so that 



4> = ln(r)S'(r), := r(ln(r)S'(r))' 



(35) 
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the potential starts out like a In potential and then smooths out, however 
becomes negative which is unphysical. 

For short distance behaviour one can again use smoothed out step func- 
tions, however numerical studies show that there are rings of large density 
near the smoothed out step. To produce a relativistic model with potential 
()33p is a substantial numerical problem: one needs to first evaluate the inte- 
gral numerically, then substitute this into the generalization of ()24p to solve 
that numerically, and finally put this into the generalizations of (j23p to get 
the pressure and density. 

Another problem is what units of length r is measured in. As galactic 
rotation curves are observed to be constant at almost all distances from the 
centre of galaxies, the problem is in this sense scale invariant and there is 
no characteristic length scale for the problem. However the models here of 
galactic rotation do have a characteristic length, usually when r = 1 where 
ln(r) changes sign, so the question arises 'what is r = 1 in meters?'. An 
exception is the model (jSSp where there is no characteristic length associated 
with a sign change of the potential, but there is the characteristic length of 
the step at Tq. There seems to be two possible answers to the problem, 
the first is that the characteristic length is the planck length, the second is 
that the characteristic length involves some new physical constant, although 
there is as yet no way of knowing how to measure it. 

6 Conclusion. 

That galactic rotation is described by a logarithmic potential, as opposed 
to a reciprocal potential — 1/r, is sufficient to show that dark matter does 
not exist. 
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